Bioactive metabolites of Asparagopsis stabilized in canola oil completely suppress methane emissions in beef cattle fed a feedlot diet.

Asparagopsis taxiformis (Asparagopsis) has been shown to be highly efficacious at inhibiting the production of methane (CH4) in ruminants. To date, Asparagopsis has been primarily produced as a dietary supplement by freeze-drying to retain the volatile bioactive compound bromoform (CHBr3) in the product. Steeping of Asparagopsis bioactive compounds into a vegetable oil carrier (Asp-Oil) is an alternative method of stabilizing Asparagopsis as a ruminant feed additive. A dose-response experimental design used 3 Asp-Oil-canola oil blends, low, medium, and high Asp-Oil which provided 17, 34, and 51 mg Asparagopsis derived CHBr3/kg dry matter intake (DMI), respectively (in addition to a zero CHBr3 canola oil control), in a tempered-barley based feedlot finisher diet, fed for 59 d to 20 Angus heifers (five replicates per treatment). On four occasions, live weight was measured and CH4 emissions were quantified in respiration chambers, and blood, rumen fluid, and fecal samples were collected. At the end of the experiment, all animals were slaughtered, with carcasses graded, and samples of meat and edible offal collected for testing of consumer sensory qualities and residues of CHBr3, bromide, and iodide. All Asp-Oil treatments reduced CH4 yield (g CH4/kg DMI, P = 0.008) from control levels, with the low, medium, and high Asp-Oil achieving 64%, 98%, and 99% reduction, respectively. Dissolved hydrogen increased linearly with increasing Asp-Oil inclusion, by more than 17-fold in the high Asp-Oil group (P = 0.017). There was no effect of Asp-Oil treatment on rumen temperature, pH, reduction potential, volatile fatty acid and ammonia production, rumen pathology, and histopathology (P > 0.10). There were no differences in animal production and carcass parameters (P > 0.10). There was no detectable CHBr3 in feces or any carcass samples (P > 0.10), and iodide and bromide residues in kidneys were at levels unlikely to lead to consumers exceeding recommended maximum intakes. Overall, Asp-Oil was found to be safe for animals and consumers of meat, and effective at reducing CH4 emissions and yield by up to 99% within the range of inclusion levels tested.

Red seaweed, Asparagopsis taxiformis (Asparagopsis), has been shown to be highly effective at inhibiting the production of methane (CH4) in ruminants. An alternative to feeding whole, freeze-dried Asparagopsis is steeping the biomass in vegetable oil to stabilize the bioactive compounds (Asp-Oil) and feeding Asp-Oil to ruminants as a component of their dietary intake. This experiment measured the CH4 reduction potential and safety of Asp-Oil in a trial with 20 Angus heifers, fed iso-fat feedlot diets containing one of the three levels of Asp-Oil, or a control oil. Compared to the control, bromoform inclusion levels of 17, 34, and 51 mg/kg of dry matter (DM; low, medium, high) reduced CH4 yield (g CH4/kg DM intake) by 64%, 98%, and 99%, respectively. There were no effects on animal production or carcass characteristics. There were no impacts on animal health, welfare, or rumen function. Carcasses were safe for human consumption, and there was no bromoform detected in any carcass samples. Overall, Asp-Oil was found to effectively reduce CH4 emissions and is safe for animals and consumers of meat and edible offal.

Pleomorphic adenoma gene1 in reproduction and implication for embryonic survival in cattle: a review.

The pleomorphic adenoma gene1 (PLAG1) encodes a DNA-binding, C2H2 zinc-finger protein which acts as a transcription factor that regulates the expression of diverse genes across different organs and tissues; hence, the name pleomorphic. Rearrangements of the PLAG1 gene, and/or overexpression, are associated with benign tumors and cancers in a variety of tissues. This is best described for pleomorphic adenoma of the salivary glands in humans. The most notable expression of PLAG1 occurs during embryonic and fetal development, with lesser expression after birth. Evidence has accumulated of a role for PLAG1 protein in normal early embryonic development and placentation in mammals. PLAG1 protein influences the expression of the ike growth factor 2 (IGF2) gene and production of IGF2 protein. IGF2 is an important mitogen in ovarian follicles/oocytes, embryos, and fetuses. The PLAG1-IGF2 axis, therefore, provides one pathway whereby PLAG1 protein can influence embryonic survival and pregnancy. PLAG1 also influences over 1,000 other genes in embryos including those associated with ribosomal assembly and proteins. Brahman (Bos indicus) heifers homozygous for the PLAG1 variant, rs109815800 (G > T), show greater fertility than contemporary heifers with either one, or no copy, of the variant. Greater fertility in heifers homozygous for rs109815800 could be the result of early puberty and/or greater embryonic survival. The present review first looks at the broader roles of the PLAG1 gene and PLAG1 protein and then focuses on the emerging role of PLAG1/PLAG1 in embryonic development and pregnancy. A deeper understanding of factors which influence embryonic development is required for the next transformational increase in embryonic survival and successful pregnancy for both in vivo and in vitro derived embryos in cattle.

The pleomorphic adenoma gene1 (PLAG1) produces PLAG1 protein which, by binding to specific regions on DNA, influences the activity of other genes that regulate many body functions. One gene is insulin-like growth factor 2 (IGF2) which controls cell metabolism and growth. The PLAG1 gene is particularly active during embryonic and fetal growth, and through IGF2 determines stature later in life. IGF2 protein is also very important in early embryonic development. This review explores the hypothesis that PLAG1 is an important determinant of embryonic survival and the establishment of pregnancy in mammals.

Genomic regions and biological pathways associated with sex-limited reproductive traits in bovine species.

Many animal species exhibit sex-limited traits, where certain phenotypes are exclusively expressed in one sex. Yet, the genomic regions that contribute to these sex-limited traits in males and females remain a subject of debate. Reproductive traits are ideal phenotypes to study sexual differences since they are mostly expressed in a sex-limited way. Therefore, this study aims to use local correlation analyses to identify genomic regions and biological pathways significantly associated with male and female sex-limited traits in two distinct cattle breeds (Brahman (BB) and Tropical Composite (TC)). We used the Correlation Scan method to perform local correlation analysis on 42 trait pairs consisting of six female and seven male reproductive traits recorded on ~1000 animals for each sex in each breed. To pinpoint a specific region associated these sex-limited reproductive traits, we investigated the genomic region(s) consistently identified as significant across the 42 trait pairs in each breed. The genes found in the identified regions were subjected to Quantitative Trait Loci (QTL) colocalization, QTL enrichment analyses, and functional analyses to gain biological insight into sexual differences. We found that the genomic regions associated with the sex-limited reproductive phenotypes are widely distributed across all the chromosomes. However, no single region across the genome was associated all the 42 reproductive trait pairs in the two breeds. Nevertheless, we found a region on the X-chromosome to be most significant for 80-90% (BB; 33 and TC; 38) of the total 42 trait pairs. A considerable number of the genes in this region were regulatory genes. By considering only genomic regions that were significant for at least 50% of the 42 trait pairs, we observed more regions spread across the autosomes and the X-chromosome. All genomic regions identified were highly enriched for trait-specific QTL linked to sex-limited traits (percentage normal sperm, metabolic weight, average daily gain, carcass weight, age at puberty, etc.). The gene list created from these identified regions were enriched for biological pathways that contribute to the observed differences between sexes. Our results demonstrate that genomic regions associated with male and female sex-limited reproductive traits are distributed across the genome. Yet, chromosome X seems to exert a relatively larger effect on the phenotypic variation observed between the sexes.

Multi-breed genomic predictions and functional variants for fertility of tropical bulls.

Worldwide, most beef breeding herds are naturally mated. As such, the ability to identify and select fertile bulls is critically important for both productivity and genetic improvement. Here, we collected ten fertility-related phenotypes for 6,063 bulls from six tropically adapted breeds. Phenotypes were comprised of four bull conformation traits and six traits directly related to the quality of the bull's semen. We also generated high-density DNA genotypes for all the animals. In total, 680,758 single nucleotide polymorphism (SNP) genotypes were analyzed. The genomic correlation of the same trait observed in different breeds was positive for scrotal circumference and sheath score on most breed comparisons, but close to zero for the percentage of normal sperm, suggesting a divergent genetic background for this trait. We confirmed the importance of a breed being present in the reference population to the generation of accurate genomic estimated breeding values (GEBV) in an across-breed validation scenario. Average GEBV accuracies varied from 0.19 to 0.44 when the breed was not included in the reference population. The range improved to 0.28 to 0.59 when the breed was in the reference population. Variants associated with the gene HDAC4, six genes from the spermatogenesis-associated (SPATA) family of proteins, and 29 transcription factors were identified as candidate genes. Collectively these results enable very early in-life selection for bull fertility traits, supporting genetic improvement strategies currently taking place within tropical beef production systems. This study also improves our understanding of the molecular basis of male fertility in mammals.

Correlation scan: identifying genomic regions that affect genetic correlations applied to fertility traits.

Although the genetic correlations between complex traits have been estimated for more than a century, only recently we have started to map and understand the precise localization of the genomic region(s) that underpin these correlations. Reproductive traits are often genetically correlated. Yet, we don't fully understand the complexities, synergism, or trade-offs between male and female fertility. In this study, we used reproductive traits in two cattle populations (Brahman; BB, Tropical Composite; TC) to develop a novel framework termed correlation scan (CS). This framework was used to identify local regions associated with the genetic correlations between male and female fertility traits. Animals were genotyped with bovine high-density single nucleotide polymorphisms (SNPs) chip assay. The data used consisted of ~1000 individual records measured through frequent ovarian scanning for age at first corpus luteum (AGECL) and a laboratory assay for serum levels of insulin growth hormone (IGF1 measured in bulls, IGF1b, or cows, IGF1c). The methodology developed herein used correlations of 500-SNP effects in a 100-SNPs sliding window in each chromosome to identify local genomic regions that either drive or antagonize the genetic correlations between traits. We used Fisher's Z-statistics through a permutation method to confirm which regions of the genome harboured significant correlations. About 30% of the total genomic regions were identified as driving and antagonizing genetic correlations between male and female fertility traits in the two populations. These regions confirmed the polygenic nature of the traits being studied and pointed to genes of interest. For BB, the most important chromosome in terms of local regions is often located on bovine chromosome (BTA) 14. However, the important regions are spread across few different BTA's in TC. Quantitative trait loci (QTLs) and functional enrichment analysis revealed many significant windows co-localized with known QTLs related to milk production and fertility traits, especially puberty. In general, the enriched reproductive QTLs driving the genetic correlations between male and female fertility are the same for both cattle populations, while the antagonizing regions were population specific. Moreover, most of the antagonizing regions were mapped to chromosome X. These results suggest regions of chromosome X for further investigation into the trade-offs between male and female fertility. We compared the CS with two other recently proposed methods that map local genomic correlations. Some genomic regions were significant across methods. Yet, many significant regions identified with the CS were overlooked by other methods.

Utility of multi-omics data to inform genomic prediction of heifer fertility traits.

Biologically informed single nucleotide polymorphisms (SNPs) impact genomic prediction accuracy of the target traits. Our previous genomics, proteomics, and transcriptomics work identified candidate genes related to puberty and fertility in Brahman heifers. We aimed to test this biological information for capturing heritability and predicting heifer fertility traits in another breed i.e., Tropical Composite. The SNP from the identified genes including 10 kilobases (kb) region on either side were selected as biologically informed SNP set. The SNP from the rest of the Bos taurus genes including 10-kb region on either side were selected as biologically uninformed SNP set. Bovine high-density (HD) complete SNP set (628,323 SNP) was used as a control. Two populations-Tropical Composites (N = 1331) and Brahman (N = 2310)-had records for three traits: pregnancy after first mating season (PREG1, binary), first conception score (FCS, score 1 to 3), and rebreeding score (REB, score 1 to 3.5). Using the best linear unbiased prediction method, effectiveness of each SNP set to predict the traits was tested in two scenarios: a 5-fold cross-validation within Tropical Composites using biological information from Brahman studies, and application of prediction equations from one breed to the other. The accuracy of prediction was calculated as the correlation between genomic estimated breeding values and adjusted phenotypes. Results show that biologically informed SNP set estimated heritabilities not significantly better than the control HD complete SNP set in Tropical Composites; however, it captured all the observed genetic variance in PREG1 and FCS when modeled together with the biologically uninformed SNP set. In 5-fold cross-validation within Tropical Composites, the biologically informed SNP set performed marginally better (statistically insignificant) in terms of prediction accuracies (PREG1: 0.20, FCS: 0.13, and REB: 0.12) as compared to HD complete SNP set (PREG1: 0.17, FCS: 0.10, and REB: 0.11), and biologically uninformed SNP set (PREG1: 0.16, FCS: 0.10, and REB: 0.11). Across-breed use of prediction equations still remained a challenge: accuracies by all SNP sets dropped to around zero for all traits. The performance of biologically informed SNP was not significantly better than other sets in Tropical Composites. However, results indicate that biological information obtained from Brahman was successful to predict the fertility traits in Tropical Composite population.

Prior biological information can be helpful in the genomic prediction of the traits. Previous multi-omics studies by our group identified genes relevant to puberty and fertility in Brahman cattle, a beef breed in northern Australia. We used this biological information in the genomic prediction of three heifer fertility traits, measured in another beef cattle breed: Tropical Composites. The three traits were: pregnancy status after the first mating season (PREG1), first conception score (FCS), and rebreeding score (REB). To test if prior biological information could capture genetic variation in the traits and improve genomic predictions, we compared the results obtained using three subsets of genetic information (i.e., subsets of DNA variants). The first subset contained only variants deemed biologically relevant (as per previous multi-omics studies). The second subset contained only variants considered biologically irrelevant. The third subset had all the variants contained in the commercial DNA assay known as the bovine high-density chip, intended as a practical control. The results indicate that multi-omics data was informative across breed scenario and can be useful in informing genomic predictions of traits of interest.

Weak genotype x environment interaction suggests that measuring scrotal circumference at 12 and 18 mo of age is helpful to select precocious Brahman cattle.

The aim of this study was to evaluate the genotype x environment interaction (GxE) for scrotal circumference (SC) measured at different ages using pedigree-based (A-1) and pedigree and genomic-based (H-1) relationship matrices. Data from 1,515 Brahman bulls, from the Cooperative Research Centre for Beef Genetic Technologies (Beef CRC) experimental dataset, were used in this study. SC was adjusted to age and body weight measured at 6 mo (SC6), 12 mo (SC12), 18 mo (SC18), and 24 mo of age (SC24). Body weight (BW) measured at 6 mo (BW6), 12 mo (BW12), 18 mo (BW18), and 24 mo of age (BW24) was used as criteria to describe the environment for SC in each age. All the animals measured were genotyped using medium-density SNP chips ("50k" or "70k" SNP) and their genotype were imputed using a reference panel with 729,068 SNP. The environment gradient (EG) was obtained by standardizing the solutions of the contemporary groups obtained by Animal Model with BW as the dependent variable. Then, the reaction norms (RN) were determined through a Random Regression Model. The breeding values (EBV) were estimated using either A-1 or H-1. The rank correlation was obtained using Spearman's correlation among the EBV estimated for the traits in analysis. For SC6 and SC24, higher estimates of heritability (h²) were obtained using A-1, when compared with those observed with H-1. In those ages, the improvement of the environment decreases the h² coefficient. On the other hand, the h² for SC12 and SC18 increased as the environment became more favorable, regardless of the matrix used. The RN for SC6 and SC24 estimated using A-1 and H-1 showed a decrease of variance from the worst to the best environment, an indication of existence of GxE. On the other hand, for SC12 and SC18, there were no significant differences between the EBV estimated in the lower and in the higher environments, regardless of the kinship matrix used, suggesting absence of GxE on those ages. Spearman's correlation among EBV estimated using A-1 and H-1 in different EG was practically equal to unity for all traits evaluated. In our study, there was weak evidence of GxE effect on SC in ages suitable for selection for sexual precocity. So, the absence of GxE at 12 and 18 mo means that these ages are advantageous for measuring SC to selection for sexual precocity. The advantage is that no changes in classification were observed when the sires were evaluated in different environments.

Beef production systems rely on efficient cow-calf operations, that is, farms where the cow herd has a high level of fertility and pregnancies are common. Bull fertility also plays an important role in terms of pregnancy rates. To increase herd fertility, cattle breeders and genetic selection programs use some indicator traits that are related to fertility. A common indicator trait used is scrotal circumference (SC), which like most reproduction traits are influenced by the animal's genetics and its environment. For some traits, when the environment has a large effect and it interacts with the genetics of the animals, selection might be less successful. Therefore, it is important to investigate genotype by environment interactions and their effect on reproduction traits used for selection. In this study, SC was measured at four different ages in Brahman cattle. We found weak evidence of genotype by environment effect on SC measured at 12 and 18 mo. In short, SC measured at these ages can be a good indicator of sexual precocity. No changes in sire rankings were observed when SC was measured at those ages, meaning that selecting the best sire is more straightforward than if the environment was playing a bigger role.

Alanine-specific appetite in slow growing chickens is associated with impaired glucose transport and TCA cycle.

BACKGROUND: The rate of protein accretion and growth affect amino acid requirements in young animals. Differences in amino acid metabolism contribute to individual variations in growth rate. This study aimed at determining how amino acid needs may change with growth rates in broiler chickens. Experiment 1 consisted of testing amino acid choices in two chicken groups with extreme growth rates (the slowest -SG- or fastest -FG- growing birds in a flock). Essential (EAA) (methionine, lysine and threonine) or non-essential (NEAA) (alanine, aspartic acid and asparagine) amino acids were added to a standard control feed (13.2 MJ/kg; 21.6% crude protein). The chickens were offered simultaneous access to the control feed and a feed supplemented with one of the two amino acid mixes added at 73% above standard dietary levels. Experiment 2 consisted of the selection of the bottom 5 SG and top 5 FG chickens from a flock of 580 to study differences in amino acid metabolism using the proventriculus representing gut sensing mechanism. In this experiment, transcriptomic, proteomic, and genomic analyses were used to compare the two groups of chickens. RESULTS: SG preferred NEAA, while they rejected EAA supplemented feeds (P < 0.05). However, FG rejected NEAA (P < 0.05), and they were indifferent to EAA supplemented feed (P > 0.05). Transcriptomic and proteomic analyses identified 909 differentially expressed genes and 146 differentially abundant proteins associated with differences in growth rate (P < 0.05). The integration of gene expression and protein abundance patterns showed the downregulation of sensing and transport of alanine and glucose associated with increased alanine catabolism to pyruvate in SG chickens. CONCLUSION: Dietary preferences for NEAA in the SG group are associated with a potential cytosolic depletion of alanine following an upregulation of the catabolism into TCA cycle intermediates.

Chromatin accessibility and regulatory vocabulary across indicine cattle tissues.

BACKGROUND: Spatiotemporal changes in the chromatin accessibility landscape are essential to cell differentiation, development, health, and disease. The quest of identifying regulatory elements in open chromatin regions across different tissues and developmental stages is led by large international collaborative efforts mostly focusing on model organisms, such as ENCODE. Recently, the Functional Annotation of Animal Genomes (FAANG) has been established to unravel the regulatory elements in non-model organisms, including cattle. Now, we can transition from prediction to validation by experimentally identifying the regulatory elements in tropical indicine cattle. The identification of regulatory elements, their annotation and comparison with the taurine counterpart, holds high promise to link regulatory regions to adaptability traits and improve animal productivity and welfare. RESULTS: We generate open chromatin profiles for liver, muscle, and hypothalamus of indicine cattle through ATAC-seq. Using robust methods for motif discovery, motif enrichment and transcription factor binding sites, we identify potential master regulators of the epigenomic profile in these three tissues, namely HNF4, MEF2, and SOX factors, respectively. Integration with transcriptomic data allows us to confirm some of their target genes. Finally, by comparing our results with Bos taurus data we identify potential indicine-specific open chromatin regions and overlaps with indicine selective sweeps. CONCLUSIONS: Our findings provide insights into the identification and analysis of regulatory elements in non-model organisms, the evolution of regulatory elements within two cattle subspecies as well as having an immediate impact on the animal genetics community in particular for a relevant productive species such as tropical cattle.

Meta-Analysis of Heifer Traits Identified Reproductive Pathways in Bos indicus Cattle.

Fertility traits measured early in life define the reproductive potential of heifers. Knowledge of genetics and biology can help devise genomic selection methods to improve heifer fertility. In this study, we used ~2400 Brahman cattle to perform GWAS and multi-trait meta-analysis to determine genomic regions associated with heifer fertility. Heifer traits measured were pregnancy at first mating opportunity (PREG1, a binary trait), first conception score (FCS, score 1 to 3) and rebreeding score (REB, score 1 to 3.5). The heritability estimates were 0.17 (0.03) for PREG1, 0.11 (0.05) for FCS and 0.28 (0.05) for REB. The three traits were highly genetically correlated (0.75-0.83) as expected. Meta-analysis was performed using SNP effects estimated for each of the three traits, adjusted for standard error. We identified 1359 significant SNPs (p-value < 9.9 × 10-6 at FDR < 0.0001) in the multi-trait meta-analysis. Genomic regions of 0.5 Mb around each significant SNP from the meta-analysis were annotated to create a list of 2560 positional candidate genes. The most significant SNP was in the vicinity of a genomic region on chromosome 8, encompassing the genes SLC44A1, FSD1L, FKTN, TAL2 and TMEM38B. The genomic region in humans that contains homologs of these genes is associated with age at puberty in girls. Top significant SNPs pointed to additional fertility-related genes, again within a 0.5 Mb region, including ESR2, ITPR1, GNG2, RGS9BP, ANKRD27, TDRD12, GRM1, MTHFD1, PTGDR and NTNG1. Functional pathway enrichment analysis resulted in many positional candidate genes relating to known fertility pathways, including GnRH signaling, estrogen signaling, progesterone mediated oocyte maturation, cAMP signaling, calcium signaling, glutamatergic signaling, focal adhesion, PI3K-AKT signaling and ovarian steroidogenesis pathway. The comparison of results from this study with previous transcriptomics and proteomics studies on puberty of the same cattle breed (Brahman) but in a different population identified 392 genes in common from which some genes-BRAF, GABRA2, GABR1B, GAD1, FSHR, CNGA3, PDE10A, SNAP25, ESR2, GRIA2, ORAI1, EGFR, CHRNA5, VDAC2, ACVR2B, ORAI3, CYP11A1, GRIN2A, ATP2B3, CAMK2A, PLA2G, CAMK2D and MAPK3-are also part of the above-mentioned pathways. The biological functions of the positional candidate genes and their annotation to known pathways allowed integrating the results into a bigger picture of molecular mechanisms related to puberty in the hypothalamus-pituitary-ovarian axis. A reasonable number of genes, common between previous puberty studies and this study on early reproductive traits, corroborates the proposed molecular mechanisms. This study identified the polymorphism associated with early reproductive traits, and candidate genes that provided a visualization of the proposed mechanisms, coordinating the hypothalamic, pituitary, and ovarian functions for reproductive performance in Brahman cattle.

Network Analyses Predict Small RNAs That Might Modulate Gene Expression in the Testis and Epididymis of Bos indicus Bulls.

Spermatogenesis relies on complex molecular mechanisms, essential for the genesis and differentiation of the male gamete. Germ cell differentiation starts at the testicular parenchyma and finishes in the epididymis, which has three main regions: head, body, and tail. RNA-sequencing data of the testicular parenchyma (TP), head epididymis (HE), and tail epididymis (TE) from four bulls (three biopsies per bull: 12 samples) were subjected to differential expression analyses, functional enrichment analyses, and co-expression analyses. The aim was to investigate the co-expression and infer possible regulatory roles for transcripts involved in the spermatogenesis of Bos indicus bulls. Across the three pairwise comparisons, 3,826 differentially expressed (DE) transcripts were identified, of which 384 are small RNAs. Functional enrichment analysis pointed to gene ontology (GO) terms related to ion channel activity, detoxification of copper, neuroactive receptors, and spermatogenesis. Using the regulatory impact factor (RIF) algorithm, we detected 70 DE small RNAs likely to regulate the DE transcripts considering all pairwise comparisons among tissues. The pattern of small RNA co-expression suggested that these elements are involved in spermatogenesis regulation. The 3,826 DE transcripts (mRNAs and small RNAs) were further subjected to co-expression analyses using the partial correlation and information theory (PCIT) algorithm for network prediction. Significant correlations underpinned the co-expression network, which had 2,216 transcripts connected by 158,807 predicted interactions. The larger network cluster was enriched for male gamete generation and had 15 miRNAs with significant RIF. The miRNA bta-mir-2886 showed the highest number of connections (601) and was predicted to down-regulate ELOVL3, FEZF2, and HOXA13 (negative co-expression correlations and confirmed with TargetScan). In short, we suggest that bta-mir-2886 and other small RNAs might modulate gene expression in the testis and epididymis, in Bos indicus cattle.

Genome-Wide Co-Expression Distributions as a Metric to Prioritize Genes of Functional Importance.

Genome-wide gene expression analysis are routinely used to gain a systems-level understanding of complex processes, including network connectivity. Network connectivity tends to be built on a small subset of extremely high co-expression signals that are deemed significant, but this overlooks the vast majority of pairwise signals. Here, we developed a computational pipeline to assign to every gene its pair-wise genome-wide co-expression distribution to one of 8 template distributions shapes varying between unimodal, bimodal, skewed, or symmetrical, representing different proportions of positive and negative correlations. We then used a hypergeometric test to determine if specific genes (regulators versus non-regulators) and properties (differentially expressed or not) are associated with a particular distribution shape. We applied our methodology to five publicly available RNA sequencing (RNA-seq) datasets from four organisms in different physiological conditions and tissues. Our results suggest that genes can be assigned consistently to pre-defined distribution shapes, regarding the enrichment of differential expression and regulatory genes, in situations involving contrasting phenotypes, time-series, or physiological baseline data. There is indeed a striking additional biological signal present in the genome-wide distribution of co-expression values which would be overlooked by currently adopted approaches. Our method can be applied to extract further information from transcriptomic data and help uncover the molecular mechanisms involved in the regulation of complex biological process and phenotypes.

X chromosome variants are associated with male fertility traits in two bovine populations.

BACKGROUND: Twenty-five phenotypes were measured as indicators of bull fertility (1099 Brahman and 1719 Tropical Composite bulls). Measurements included sperm morphology, scrotal circumference, and sperm chromatin phenotypes such as DNA fragmentation and protamine deficiency. We estimated the heritability of these phenotypes and carried out genome-wide association studies (GWAS) within breed, using the bovine high-density chip, to detect quantitative trait loci (QTL). RESULTS: Our analyses suggested that both sperm DNA fragmentation and sperm protamine deficiency are heritable (h2 from 0.10 to 0.22). To confirm these first estimates of heritability, further studies on sperm chromatin traits, with larger datasets are necessary. Our GWAS identified 12 QTL for bull fertility traits, based on at least five polymorphisms (P < 10-8) for each QTL. Five QTL were identified in Brahman and another seven in Tropical Composite bulls. Most of the significant polymorphisms detected in both breeds and nine of the 12 QTL were on chromosome X. The QTL were breed-specific, but for some traits, a closer inspection of the GWAS results revealed suggestive single nucleotide polymorphism (SNP) associations (P < 10-7) in both breeds. For example, the QTL for inhibin level in Braham could be relevant to Tropical Composites too (many polymorphisms reached P < 10-7 in the same region). The QTL for sperm midpiece morphological abnormalities on chromosome X (QTL peak at 4.92 Mb, P < 10-17) is an example of a breed-specific QTL, supported by 143 significant SNPs (P < 10-8) in Brahman, but absent in Tropical Composites. Our GWAS results add evidence to the mammalian specialization of the X chromosome, which during evolution has accumulated genes linked to spermatogenesis. Some of the polymorphisms on chromosome X were associated to more than one genetically correlated trait (correlations ranged from 0.33 to 0.51). Correlations and shared polymorphism associations support the hypothesis that these phenotypes share the same underlying cause, i.e. defective spermatogenesis. CONCLUSIONS: Genetic improvement for bull fertility is possible through genomic selection, which is likely more accurate if the QTL on chromosome X are considered in the predictions. Polymorphisms associated with male fertility accumulate on this chromosome in cattle, as in humans and mice, suggesting its specialization.

Dynamics of Gene Co-expression Networks in Time-Series Data: A Case Study in Drosophila melanogaster Embryogenesis.

Co-expression networks tightly coordinate the spatiotemporal patterns of gene expression unfolding during development. Due to the dynamic nature of developmental processes simply overlaying gene expression patterns onto static representations of co-expression networks may be misleading. Here, we aim to formally quantitate topological changes of co-expression networks during embryonic development using a publicly available Drosophila melanogaster transcriptome data set comprising 14 time points. We deployed a network approach which inferred 10 discrete co-expression networks by smoothly sliding along from early to late development using 5 consecutive time points per window. Such an approach allows changing network structure, including the presence of hubs, modules and other topological parameters to be quantitated. To explore the dynamic aspects of gene expression captured by our approach, we focused on regulator genes with apparent influence over particular aspects of development. Those key regulators were selected using a differential network algorithm to contrast the first 7 (early) with the last 7 (late) developmental time points. This assigns high scores to genes whose connectivity to abundant differentially expressed target genes has changed dramatically between states. We have produced a list of key regulators - some increasing (e.g., Tusp, slbo, Sidpn, DCAF12, and chinmo) and some decreasing (Rfx, bap, Hmx, Awh, and mld) connectivity during development - which reflects their role in different stages of embryogenesis. The networks we have constructed can be explored and interpreted within Cytoscape software and provide a new systems biology approach for the Drosophila research community to better visualize and interpret developmental regulation of gene expression.

Gene regulation could be attributed to TCF3 and other key transcription factors in the muscle of pubertal heifers.

Puberty is a whole-body event, driven by the hypothalamic integration of peripheral signals such as leptin or IGF-1. In the process of puberty, reproductive development is simultaneous to growth, including muscle growth. To enhance our understanding of muscle function related to puberty, we performed transcriptome analyses of muscle samples from six pre- and six post-pubertal Brahman heifers (Bos indicus). Our aims were to perform differential expression analyses and co-expression analyses to derive a regulatory gene network associate with puberty. As a result, we identified 431 differentially expressed (DEx) transcripts (genes and non-coding RNAs) when comparing pre- to post-pubertal average gene expression. The DEx transcripts were compared with all expressed transcripts in our samples (over 14,000 transcripts) for functional enrichment analyses. The DEx transcripts were associated with "extracellular region," "inflammatory response" and "hormone activity" (adjusted p < .05). Inflammatory response for muscle regeneration is a necessary aspect of muscle growth, which is accelerated during puberty. The term "hormone activity" may signal genes that respond to progesterone signalling in the muscle, as the presence of this hormone is an important difference between pre- and post-pubertal heifers in our experimental design. The DEx transcript with the highest average expression difference was a mitochondrial gene, ENSBTAG00000043574 that might be another important link between energy metabolism and puberty. In the derived co-expression gene network, we identified six hub genes: CDC5L, MYC, TCF3, RUNX2, ATF2 and CREB1. In the same network, 48 key regulators of DEx transcripts were identified, using a regulatory impact factor metric. The hub gene TCF3 was also a key regulator. The majority of the key regulators (22 genes) are members of the zinc finger family, which has been implicated in bovine puberty in other tissues. In conclusion, we described how puberty may affect muscle gene expression in cattle.

Polymorphisms in major histocompatibility complex genes and its associations with milk quality in Murrah buffaloes.

Animal breeding programs have used molecular genetic tools as an auxiliary method to identify and select animals with superior genetic merit for milk production and milk quality traits as well as disease resistance. Genes of the major histocompatibility complex (MHC) are important molecular markers for disease resistance that could be applied for genetic selection. The aim of this study was to identify single nucleotide polymorphisms (SNPs) and haplotypes in DRB2, DRB3, DMA, and DMB genes in Murrah breed and to analyze the association between molecular markers and milk, fat, protein and mozzarella production, fat and protein percentage, and somatic cell count. Two hundred DNA samples from Murrah buffaloes were used. The target regions of candidate genes were amplified by polymerase chain reaction (PCR) followed by sequencing and identification of polymorphisms. Allele and genotype frequencies, as well as linkage disequilibrium between SNPs, were calculated. Genotypes were used in association analyses with milk production and quality traits. Except for the DMA gene, identified as monomorphic, the other genes presented several polymorphisms. The DMB, DRB2, and DRB3 genes presented two, six, and seven SNPs, respectively. Fifty-seven haplotype blocks were constructed from 15 SNPs identified, which was used in association analyses. All the studied traits had at least one associated haplotype. In conclusion, it is suggested that the haplotypes found herein can be associated with important traits related to milk production and quality.

Proteomics Recapitulates Ovarian Proteins Relevant to Puberty and Fertility in Brahman Heifers (Bos indicus L.).

High fertility and early puberty in Bos indicus heifers are desirable and genetically correlated traits in beef production. The hypothalamus-pituitary-ovarian (HPO) axis synthesizes steroid hormones, which contribute to the shift from the pre-pubertal state into the post-pubertal state and influence subsequent fertility. Understanding variations in abundance of proteins that govern steroid synthesis and ovarian signaling pathways remains crucial to understanding puberty and fertility. We used whole ovaries of six pre-pubertal and six post-pubertal Brahman heifers to conduct differential abundance analyses of protein profiles between the two physiological states. Extracted proteins were digested into peptides followed by identification and quantification with massspectrometry (MS) by sequential window acquisition of all instances of theoretical fragment ion mass spectrometry (SWATH-MS). MS and statistical analysis identified 566 significantly differentially abundant (DA) proteins (adjusted p < 0.05), which were then analyzed for gene ontology and pathway enrichment. Our data indicated an up-regulation of steroidogenic proteins contributing to progesterone synthesis at luteal phase post-puberty. Proteins related to progesterone signaling, TGF-ß, retinoic acid, extracellular matrix, cytoskeleton, and pleiotrophin signaling were DA in this study. The DA proteins probably relate to the formation and function of the corpus luteum, which is only present after ovulation, post-puberty. Some DA proteins might also be related to granulosa cells signaling, which regulates oocyte maturation or arrest in ovaries prior to ovulation. Ten DA proteins were coded by genes previously associated with reproductive traits according to the animal quantitative trait loci (QTL) database. In conclusion, the DA proteins and their pathways were related to ovarian activity in Bos indicus cattle. The genes that code for these proteins may explain some known QTLs and could be targeted in future genetic studies.

RAPID COMMUNICATION: Multi-breed validation study unraveled genomic regions associated with puberty traits segregating across tropically adapted breeds1.

An efficient strategy to improve QTL detection power is performing across-breed validation studies. Variants segregating across breeds are expected to be in high linkage disequilibrium (LD) with causal mutations affecting economically important traits. The aim of this study was to validate, in a Tropical Composite cattle (TC) population, QTL associations identified for sexual precocity traits in a Nellore and Brahman meta-analysis genome-wide association study. In total, 2,816 TC, 8,001 Nellore, and 2,210 Brahman animals were available for the analysis. For that, genomic regions significantly associated with puberty traits in the meta-analysis study were validated for the following sexual precocity traits in TC: age at first corpus luteum (AGECL), first postpartum anestrus interval (PPAI), and scrotal circumference at 18 months of age (SC). We considered validated QTL those underpinned by significant markers from the Nellore and Brahman meta-analysis (P ≤ 10-4) that were also significant for a TC trait, i.e., presenting a P-value of ≤10-3 for AGECL, PPAI, or SC. We also considered as validated QTL those regions where significant markers in the reference population were at ±250 kb from significant markers in the validation population. Using this criteria, 49 SNP were validated for AGECL, 4 for PPAI, and 14 for SC, from which 5 were in common with AGECL, totaling 62 validated SNP for these traits and 30 candidate genes surrounding them. Considering just candidate genes closest to the top SNP of each chromosome, for AGECL 8 candidate genes were identified: COL8A1, PENK, ENSBTAG00000047425, BPNT1, ADAMTS17, CCHCR1, SUFU, and ENSBTAG00000046374. For PPAI, 3 genes emerged as candidates (PCBP3, KCNK10, and MRPS5), and for SC 8 candidate genes were identified (SNORA70, TRAC, ASS1, BPNT1, LRRK1, PKHD1, PTPRM, and ENSBTAG00000045690). Several candidate regions presented here were previously associated with puberty traits in cattle. The majority of emerging candidate genes are related to biological processes involved in reproductive events, such as maintenance of gestation, and some are known to be expressed in reproductive tissues. Our results suggested that some QTL controlling early puberty seem to be segregating across cattle breeds adapted to tropical conditions.

Proteomics Reveals Profound Metabolic Changes in the Alcohol Use Disorder Brain.

Changes in brain metabolism are a hallmark of alcohol use disorder (AUD). Determining how AUD changes the brain proteome is critical for understanding the effects of alcohol consumption on biochemical processes in the brain. We used data-independent acquisition mass spectrometry proteomics to study differences in the abundance of proteins associated with AUD in prefrontal lobe and motor cortex from autopsy brain. AUD had a substantial effect on the overall brain proteome exceeding the inherent differences between brain regions. Proteins associated with glycolysis, trafficking, the cytoskeleton, and excitotoxicity were altered in abundance in AUD. We observed extensive changes in the abundance of key metabolic enzymes, consistent with a switch from glucose to acetate utilization in the AUD brain. We propose that metabolic adaptations allowing efficient acetate utilization contribute to ethanol dependence in AUD.

Multitrait meta-analysis identified genomic regions associated with sexual precocity in tropical beef cattle.

Multitrait meta-analyses are a strategy to produce more accurate genome-wide association studies, especially for complex phenotypes. We carried out a meta-analysis study for traits related to sexual precocity in tropical beef cattle (Nellore and Brahman) aiming to identify important genomic regions affecting these traits. The traits included in the analyses were age at first calving (AFC), early pregnancy (EP), age at first corpus luteum (AGECL), first postpartum anoestrus interval (PPAI), and scrotal circumference (SC). The traits AFC, EP, and SCN were measured in Nellore cattle, while AGECL, PPAI, and SCB were measured in Brahman cattle. Meta-analysis resulted in 108 significant single-nucleotide polymorphisms (SNPs), at an empirical threshold P-value of 1.39 × 10-5 (false discovery rate [FDR] < 0.05). Within 0.5 Mb of the significant SNP, candidate genes were annotated and analyzed for functional enrichment. Most of the closest genes to the SNP with higher significance in each chromosome have been associated with important roles in reproductive function. They are TSC22D2, KLF7, ARHGAP29, 7SK, MAP3K5, TLE3, WDR5, TAF3, TMEM68, PPP1R15B, NR2F2, GALR1, SUFU, and KCNU1. We did not observe any significant SNP in BTA5, BTA12, BTA17, BTA18, BTA19, BTA20, BTA22, BTA23, BTA25, and BTA28. Although the majority of significant SNPs are in BTA14, it was identified significant associations in multiple chromosomes (19 out of 29 autosomes), which is consistent with the postulation that reproductive traits are complex polygenic phenotypes. Five proposed association regions harbor the majority of the significant SNP (76%) and were distributed over four chromosomes (P < 1.39 × 10-5, FDR < 0.05): BTA2 (5.55%) from 95 to 96 Mb, BTA4 (5.55%) from 94.1 to 94.8 Mb, BTA14 (59.26%) from 24 to 25 Mb and 29 to 30 Mb, and BTA21 (5.55%) from 6.7 Mb to 11.4 Mb. These regions harbored key genes related to reproductive function. Moreover, these genes were enriched for functional groups associated with immune response, maternal-fetal tolerance, pregnancy maintenance, embryo development, fertility, and response to stress. Further studies including other breeds and precocity traits could confirm the importance of these regions and identify new candidate regions for sexual precocity in beef cattle.